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Abstract 

We study the interaction of mutually non-interacting Klein-Gordon particles with lo- 
, calized sources on stochastically complete Riemannian surfaces. This asymptotically free 

theory requires regularization and coupling constant renormalization. Renormalization is 
. performed non-perturbatively using the orthofermion algebra technique and the princi- 

. pal operator $ is found. The principal operator is then used to obtain the bound state 

ff^ . spectrum, in terms of binding energies to single Dirac-delta function centers. The heat 



kernel method allows us to generalize this procedure to compact and Cartan-Hadamard 
type Riemannian manifolds. We make use of upper and lower bounds on the heat kernel 



■ to constrain the ground state energy from below thus confirming that our neglect of pair 

creation is justified for certain ranges of parameters in the problem. 

■ 

X 

a 1 Introduction 

Quantum field theories only make sense after a regularization and renormalization procedure 
is applied to them. However, the complexity of quantum field theories, QCD in particular, 
combined with the lack of exact solutions in most cases obscures the concept of renormalization. 
It is extremely useful for understanding the effects of renormalization in QCD to investigate 
theories that are simple yet capture some of the important features of QCD such as asymptotic 
freedom and (classical) scale invariance. 

From the above point of view, one would wish to understand renormalization in detail in 
the quantum mechanical context where the dynamics are much simpler. This has been done 
for the interaction of non-relativistic bosons with Dirac-delta function type potentials in dimen- 
sions 2 and 3 [1] where renormalization is required as well as for the relativistic Lee model [2]. 
These models all share the property of being asymtotically free. Furthermore, the interaction of 
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non-relativistic bosons with Dirac-delta potentials in 2 dimensions, as the coupHng constant is 
dimensionless, is a scale invariant theory classically. 

Berezin and Fadeev [3] were the first to investigate the interaction of non-relativistic bosons 
with Dirac-delta function potentials, this was then taken up in J. Hoppe's doctoral dissertation 
|1] in detail. Regularization and renormalization were performed by a scattering analysis O 
El [71 El [9l Hn] and by thinking of the Dirac-delta function as the limit of different potentials 
fn\ m\ . renormalization group properties of the problem were analyzed [131 HH HSj and it 
was put in a broader mathematical framework in the theory of self-adjoint extensions of non- 
interacting hamiltonians [16] . The theory we investigate in this paper is a relativistic version of 
the interaction of non-relativistic bosons with Dirac-delta function potentials. We hope to get a 
bit closer to quantum field theory by allowing for a relativistic dispersion relation for the bosons 
even though we continue to neglect particle production due to the interaction with the source 
as the hamiltonian under scrutiny conserves the number of bosons. 

We use the principal operator technique introduced in an unpublished paper by Rajeev [T7| to 
renormalize the theory nonperturbatively and to obtain its bound state spectrum. Then, along 
the lines of [1] we generalize the renormalization procedure to compact and Cartan-Hadamard 
type Riemannian manifolds in the hope of gaining some insight into how renormalization works 
and how it modifies the theory, if at all. We reach the conclusion that although the renormal- 
ization procedure that is crucial for obtaining physically meaningful results involves advanced 
mathematical techniques, it does not seem to contradict or even alter in any fundamental way 
our intuition from elementary quantum mechanics. 

The paper is organized as follows: In section [21 we write down the cut-off hamiltonian 
and explain how to obtain its spectrum through the eigenstates of the principal operator with 
vanishing eigenvalue. The principal operator is then renormalized in a way that is valid on 
Riemannian manifolds using the heat kernel. Subsection 13.11 in section [3] shows the equivalence 
between the renormalization done using the heat kernel and the usual method in field theory 
of placing a momentum space cut-off. The momentum space wavefunction of single particle 
fields and the corresponding configuration space "wavefunctions" are obtained in subsection 13.31 
of this section. The section concludes with an approximate calculation of bound state energies 
in the tunnelling regime in subsection 13.41 Perhaps the most important results of the paper 
are the lower bounds we place on the ground state energy in section [H on flat, compact and 
Cartan-Hadamard type Riemannian manifolds. 

2 General Formalism 

We consider relativistic Klein-Gordon particles interacting with a finite number of localized 
sources on general two dimensional, stochastically complete Riemannian manifolds {Ai,g). Our 
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starting point is the following Hamiltonian 

where 

27r)2V 

HT' = -Y.3i I dgxK,/2{a^,x)<j)^-\x) j dgyK,,2{a,,y)<p'^+\y) (1) 
i=i ^ 

where Qi stands for the coupling constant. In the case of a non-relativistic particle interacting 

with Dirac-delta function type potentials in 1-dimension the coupling constant should be chosen 

positive (with the above convention) for bound states to exist. However, in our case, the sign 

of the coupling constant has little significance and can be of either sign depending on how it 

is renormalized. Here N specifies the number of Dirac-delta function centers. The subscript e 

indicates the regularization that is required to make the Hamiltonian above well-defined. In the 

e — s> 0+ limit that we are interested in, the regularized Hamiltonian has a Dirac-delta function 

type point interaction and the interaction above written in terms of the heat kernel K^{x,y) is 

just a smeared out version of it. A detailed discussion on this will be given later. The positive 

frequency part (f)^~\x) of the bosonic field is given in terms of the creation operator indexed by 

(T which has the interpretation of momentum vector on flat space as 



E 



X 

'UJ, 



= \{a) + m' (2) 



is the Hermitian conjugate of this operator. The fa{x) G L^(A^) are the orthonormal 
complete eigenfunctions of the Laplace-Beltrami operator which is defined on local coordinates 
X = {x^, x"^) as 

a, 13=1 

where Qajj are the elements of the metric tensor matrix g and g"'^ are those of its inverse. 
Eigenfunctions of the Laplace operator have the following properties 

-VlUx) = X{a)Mx) , 

/ 4"W:(a;)/a'(x) =5^^, , 
Jm 

E/^(^)/'^(^) = '^J"H^-y) • (4) 

The sums in the first line of Eq. ([2]) and in the last relation of Eq. (jlj) should be replaced by 
integrals in the case of fiat and Cartan-Hadamard type Riemannian manifolds since the spectrum 



3 



of the Laplace operator on such manifolds is continuous. We should emphasize that whenever 
we talk about the eigenvalues of the Laplace operator we mean those of — so that they are 
nonnegative. 

This simple looking problem leads to divergences and some kind of regularization is manda- 
tory to get physically sensible results. In order to renormalize the Hamiltonian in eq. ([1]) we 
will use the orthofermion algebra technique developed by Rajeev [I7j. In this method, we intro- 
duce fictitous operators x ^i-iid named angels which are reminiscent of the hopping operators 
one encounters in solid state physics. These operators commute with the bosonic creation and 
annihilation operators and satisfy the relations given below 

Xix] = ^ij'^o , 

XiXj ; 

J2xlx^ = Yl^^ = T^l■ (5) 

i i 

It is seen from the second relation and its Hermitian conjugate that these operators generate a 
two-state system when acting on the Fock space. IIq and Ilj are projection operators onto the 
spaces with no angels and with one angel of type i respectively. The physical space will be the 
direct product of the space of angels and the bosonic one. We define a new operator which in 
matrix form is as follows 



G-EUn 



{Ho-E)Uo dgxK^i2{ai,x)(t)^ Kx)xi 



a 6t 
h d 

with the resolvent or the Green's function defined as 



(6) 



iG-EIUr.[y'l) (7) 

The projection of this Green's function matrix on to the no angel subspace can be written in 
two alternative ways: 



a = {a - bU-^b)-^ = {H - E)-^ 
a + a" 



'1 ■ --ifet^-i^^-i 



where the principal operator matrix $ is given by 

= - ba-^b^ (9) 

The first relation in eq. ([H]) can easily be verified using the definition of IIq given in eq. 
This shows that the projection of the resolvent of the new operator onto the no angel subspace 



4 



reproduces the Green's function of the original Hamiltonian. The reason we can project the 
Hamiltonian on to spaces with a fixed number of bosons is because it commutes with the number 
operator N = J (|;;^«pap. The spectrum, i.e. the poles, of the original Hamiltonian corresponds 
to zeros of the principal operator matrix. The difficult task of finding the spectrum of the original 
Hamiltonian thus boils down to finding the eigenstates of this finite dimensional matrix $ with 
eigenvalue zero. 

When we use eq. to find the explicit form of the principal operator matrix we get 

N N t 

I / ^ V V ■ 



i=l i=l a V'^'^ 



(^0 - E) 

J2 [ d,yK,;,{a„y)Y,^fliy) (10) 



In the second term, angel operators can be moved around freely. However, the creation and 
annihilation operators do not commute with the free Hamiltonian and this expression is modified 
as we move the annihilation operator to the right of all the creation operators and vice versa. 
After normal-ordering we get 

^1 f f 

i=l t,j a J -J 



1 x\Xj ^ [ A ( fxi^)(^\ 



a— ' ' / dgyK,/2{ai,y) 11 

[Ho - E + UJ^ + UJx) J ^/uJ^ 

In the angel formulation, single boson states correspond to the ground state in the bosonic 
subspace. To illustrate this we feed in a 1-boson state to eq. 

^:'ba~' (|0 > (g) b >b) oc ^:'b (|0 > (g) b >b) oc c^lxi > (g) |0 >b^ (12) 

Besides, the application of the annihilation operator thus the free Hamiltonian on such states 
gives zero i.e. 

(N \ N 

E ^il^^ > ) ® 1° E ^^1^^ > (g) «p|0 >B= (13) 
i=l / i=l 

Inspection of eq. flTTl) . in view of the above, shows that normal-ordering disentangles the effect 
of the bosonic operators and angels in the case of 1-particle states. In spite of the nonlinear 
realization of renormalization, it can be shown that the many particle states are given as direct 
products of the single particle ones just like in ordinary quantum mechanics. Therefore, knowl- 
edge of the spectrum of single particle states is equivalent to a knowledge of the full spectrum 
of the theory. This allows us to ignore the third term in eq. ffTTj) from now on. 
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We would like to write the principal operator in a form that is valid on Riemannian manifolds. 
We commence by separating the denominator using the identity below 



1 



1 



ujp{-E + Up) 



1 



1 



(14) 



[-E) [ujp {-E + ujp)_ 
followed by turning these denominators into exponentials with the help of the following relation: 

1 1 

' ^ • (15) 



'-E + 



pj JO 



^su^, [1 - e^^] ds 



This relation is valid provided the real parts of the terms in the exponentials satisfy 3? {up — E) > 
0. One should think of the energy as a complex variable since the principal operator relevant 
for scattering can be obtained from this one by an appropriate analytical continuation in the 
complex energy plane. 

N ^ 

= ^—xlxi-^Yl / dgxK,/2{ai,x)fa{x) / dgyK,/2{aj,y)f*{y) 

i=l ^' i,7 ^ 



ds e 



1-e 



sE 



-E) 



xlxj 



Another ingredient that we need is the subordination identity given below 



-sA 



2v^ 



s2/(4u)--uyl2 



du 

3/2 



U 



(16) 



(17) 



Utilizing these relations in the expression for the principal operator we end up with the following 
expression 



y^—x\xi-7rr^ I dgxK,/2iai,x) [ dgyK,/2iaj,y) [ 
i=i 9i ^V^r J J Jo 

POO 

J2e-^'^'^^Ux)f:iy) / dsse- 

Jo 



du 

3/2 



U 



7(4«) 



1-e 



-E) 



xlxj 



We make a short digression to give information on the heat kernel Ku{x, y), the minimal positive 
solution to the heat equation, 



d 

—Ku{x,y) = V]Ku{x,y) 
It has this eigenfunction expansion on compact Riemannian manifolds [HI [I9] 

ir„(x,y) = 5^e-"^M/.(x)/;(y) 

and the following properties on general Riemannian manifolds 

Ku{x,y)= Ku{y,x) symmetry property, 
\\m Ku{x,y) = 6g{x,y) initial condition, 

/ dgZ {x, z) Ku2{z,y) = Ku^+U2{x,y) reproducing property. 
Jm 
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(19) 
(20) 



(21) 



We assume that the heat kernel satisfies the stochastic completeness relation given below 

Ku{x,y)dgX = l (22) 



M 



This is an immediate consequence of the eigenfunction expansion, however it will hold even when 
an eigenfunction expansion may be lacking. We identify the eigenfunction expansion of the heat 
kernel in eq. f|T8l) and scale the s variable to put it in the form below 



POO 

Ku{x, y) I ds se~^ 
Jo 



/2[ai,x 



j dgyK,/2iaj,y) ^ 



-E) 



xlxj 



(23) 



We first use the symmetry property of the heat kernel, 

j~^9i Sv/tF^Jo Vu Jo 



ds se-'"/' 



1 — gsEv^ 

i-E) 



dgX / dgyK,/2{ai,x)Kuix,y)K,/2{y,aj)xlxj 



(24) 



then use the reproducing property and do an integration by parts to put the principal operator 
in the following form 

N 



-| -| poo POO 

tt9i v^t/^/o Jo 



u+e {(^i ) (^j ) Xi Xj 



(25) 



As we will discuss in more detail later, only the diagonal part of the principal operator leads to 
a divergence in the e ^ 0"^ limit. We separate the divergent and finite parts of the principal 
operator and shift the u integral keeping only the terms necessary to cut-off the divergence to 
write the principal operator as 



^ ' 1 1 /'OO noo 

= y ^ / ds e-^'/^ / du e'^^ e-'^'^'KJai. 

[9i Jo Je 

- —Y dse-'"'^ rfMe^^v^e-"™'K,(a„a,)xl 
V^r 7/ Jo Jo 



Xi Xi 



(26) 



We choose the couphng constant in the following way so as to regularize the principal operator 

-| -I foo poo 

— r = ^/ dse-''/^ du e'"^^ e-''"''Ku{ai,ai) (27) 

9i [fJ'i, e) V^r Jo Je 

This is a good place to make a digression and determine how the coupling constant changes with 

energy scale. In order to do this, we choose the renormalized coupling constant Qr^i in terms of 

the bare coupling constant go^i as 

1 1 



1 noo poo 

A Jo Je 



(28) 



and impose the (unphysical) renormalization condition 

1 



— = (29) 
The physically relevant case will be recovered if this renormalization scale Mj is chosen to be 

th 

the binding energy to the single i Dirac-delta center. The condition above determines the 
bare coupling constant as a function of the renormalization scale and the renormalized coupling 
constant becomes the following: 



oo 



9r,i 

The i beta function l3i{gr,i) is defined as follows: 



(30) 



Using this definition and imposing at the end the physical renormalization condition by setting 
Mi = /ij, we find the result given below for the beta function 

Pi{gr,^) = -^^l^lAtil r se-'"l^ds r e-"™'y^ e^'^v^ ir„(a„ a,) du (32) 



With the natural choice of the coupling constant as in eq. fl27|) . the experimentally measured 
binding energy E coincides with the bound state energy in the single Dirac-delta center case. 
Inserting the expression for the coupling constant in eq. (1271) in the previous one yields the 
regularized form of the full principal operator: 

N 



$ = ^ V / dse-''/^ / due-''"'' KJtti, Ui 
A ^ io Jo 

-| poo POO 

- ^ V / dse-''/^ due-''"''K^{ai,aj)e^^'xl 
V ^ , , Jo Jo 



x\Xi 



Xj 



^n^Jo Jo {u,u,f'' 

dgX(P^-\x)K^,{x,a,)e-'^''°-'^'> [ dgyK^,{a^,y)<P^^\y)x^XJ (33) 



In the above equation, for completeness we included the term we had dropped earlier. As this 
expression is finite in the limit e ^ 0^, we dropped the subscript e in $. This is explained in 
more detail in eq. fl5T]) and the paragraph following it. 

What we have done so far applies to Riemannian manifolds in general. However, we will 
give up on generality for a while for the sake of exposing the results more clearly. As will be 
seen, the general features regarding the renormalization of the coupling constant, behavior of 
the off-diagonal elements at large separations of the Dirac-delta centers and how the elements of 
the principal operator matrix vary with bound state energy hold on Riemannian manifolds, too. 
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3 Renormalization and Spectrum of the Theory 



3.1 Coupling Constant Renormalization on 

In this section, we will illustrate how the regularization and renormalization procedure valid on 
Riemannian manifolds that we descibed in the previous section is in fact equivalent to the usual 
momentum cut-off prescription we encounter in quantum field theory. Identifying the /p(x)'s 
with the usual Fourier exponents and replacing the heat kernels with the Dirac-delta functions, 
eq. (fTTl) becomes the following: 



TV 



E 

i=l 

E 



1 

|p|<A 



. :-^+^p) 



Xi Xi 



E 



|p|<A 



'xlxj 



|p|<A 7|k|<A V^^^^ {Ho - E + UJp + UJy, 



-«P XiXj 



(34) 



It is obvious that the second term in the first line of the above expression diverges for n > 2. This 
divergence must be regularized by placing a cut-off A in the upper limit of that integral. This 
removes the problem of infinities, however it raises a new difficulty that all physical observables 
depend on A. This is circumvented by letting the coupling constant depend on the cut-off 
momentum i.e. Qi = gi{A) in such a way that all physical observables are independent of it in the 
limit A — > oo. This is what in quantum field theory is commonly referred to as renormalization. 
We propose the simple form below for the coupling constant to ensure this independence. 



1 



|p|<A 



(2^ 



1 



(35) 



Let us now evaluate these integrals to get an explicit expression for the principal operator 
matrix on flat surfaces. First, we will evaluate the diagonal terms. Separating the denominators 
as in eq. f|T^ and using the fact that the integrands are rotationally invariant makes their 
calculation trivial. The result is (with no summation over repeated indices) 

-A 



lim 

A^oo 



27r 



de 



i2ny 



1 



_ujp{-iJ,i + ujp) iUp{-E + ujp)_ 



lim — 

A^oo 271 



In 



27r 



m 



E 



m 



m- Hi 
+ 0(1/ A) 



In 



VA2 + m^-E 
m — E 



(36) 



The choice eq. (1351) for the coupling constant in some sense absorbs the cut-off momentum depen- 
dence and the final result given above for the diagonal elements are, as promised, independent of 
this cut-off momentum in the limit A ^ oo. As will be shown momentarily, off-diagonal matrix 
elements are already finite so the cut-off momentum dependence has successfully been removed. 
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Interestingly enough, the elements of the principal operator matrix that are now finite do 
not depend on the coupling constant! The dimensionless coupling constant gi has been traded 

th 

for a dimensionful parameter /Xj which is the experimentally measured binding energy to the i 
Dirac-delta function center. This is the analog of dimensional transmutation in the field theory 
context. 

Let us now calculate the off-diagonal elements. We use the trick of separating the denominator 
of two terms into two denominators, followed by turning these fractions into integrals over some 
variable as in eq. (fT5|) 



l^/Tiy—tj) Jq Jo u ' 



p 



The separation vector is defined as dij = — a^. After grinding through the algebra we arrive 
at the simple result given below: 

,,.zi^^-*.H."«.-..]^^ (38) 

where dij is shorthand for the length of the separation vector. With a little effort, the diagonal 
elements of the principal operator can alternatively be obtained starting from the equation above. 
For this purpose, one may find useful the following expansion of the gamma function near its 
pole at zero. 

/•oo 1 

lim / x("-i)e-''^rfx = lim a-"r(n) ~--lna-7 + C(n) (39) 

Concentrating on states that contain only a single boson let us make some observations based 
on the results for the elements of the principal operator matrix. First of all, the divergence in 
the coupling constant is a mild, logarithmic one as one would expect from naive dimensional 
analysis since the coupling constant is dimensionless. Moreover, in the limit the minimum sepa- 
ration between the different Dirac-delta function centers is much greater than the Compton 
wavelength of the particle that is md™" ^ 1, off-diagonal elements of the principal operator 
matrix are much smaller than the diagonal elements and can be considered as a perturbation. 
In order to make this explicit, we separate the principal operator matrix into a diagonal (^o{E) 
and an off-diagonal 5^{E) part as follows: 

^{E) = ^d{E) + 5^{E) (40) 

Then we look for solutions of the form 

\vi >= > +\Svi > (41) 

with |t>° > satisfying 

^r,{E = fii)\v^>=0 (42) 
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When we do this we find that the variation in the binding energies due to the presence of other 
Dirac-delta function centers is given by the equation below 

This equation is reminiscent of the 2"*^ order non-degenerate perturbation theory result, however 
unlike that equation this one is correct nonperturbatively, too. Although the change in the 
binding energy to a single Dirac-delta function center is assumed to be a perturbation i.e. \6Ei\ <^ 
nothing prohibits the initial binding energy itself from being nonperturbatively strong that 
is m — fii m. The above equation is valid only if there is no degeneracy in the initial binding 
energies. A similar equation can be derived for the degenerate case. 

It is possible, in this regime, to evaluate the off-diagonal elements approximately by hand 
which becomes exact in the limit. One can then plug this result in eq. ( H3l) to find the change 
in the binding energy. For this purpose, we may start with eq. (138|) . however this way of 
finding the approximate form of the off-diagonal elements does not lend itself to generalization 
to Riemannian manifolds. Instead, we will use the form of the principal operator which we 
derived in eq. (l33ll that is valid on Riemannian manifolds. Here we just state the results that 
will be derived in subsection 13.41 



<i>ij{E - fii) y= ■ . 1— ^ i 1 / > ^-l/2 // „ // 1 y^^) 



and 



<l>ij{E = fii) = ^ — — (mdij) < -/ii/m < 1 (45) 

zn mdij\fii/m\ 

For example, in the double Dirac-delta center case the binding energies /xi > /i2 > are modified 
as follows: 



^ 1 /m-/ii exp (-2mdi2v/l-(/iiM) 



m m 27T ym + fii rndi2\n(^^ 

\m—fj.i 

E2 fi2 1 /m-/i2 exp (-2mdi2v/l- (^2/ 



m m 27r Vm + /i2 ^dulnl^^) 



(46) 



It should also be noted that whereas the absolute value of the off-diagonal elements decrease 
with decreasing energy, the on-diagonal elements increase with decreasing energy. This obser- 
vation will play an important role when we look for a lower bound on the ground state energy. 
It turns out that all of these properties of the principal operator matrix hold on Riemannian 
manifolds as well. 
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3.2 Renormalization on Riemannian Manifolds 



The renormalized principal operator appropriate for 1-particle states written in terms of the heat 
kernel after normal-ordering is given below 



I poo pc 



du 



-I POO POO 

V^^Jo Jo 



ai,aj)xlXj (47) 



Superficially, the s-integral in the first two terms of the above equation looks like it might be 
divergent at large values of u, however, this suspicion is groundless. One can convince oneself 
that this is the case by completing the square and shifting the s-integral. Eq. fl20|) implies that 
the asymptotic form of the heat kernel at late timeil] is given by 

lim Ku{x,y) oc e""^™" (48) 

where Xmin is the minimum eigenvalue of the Laplacian on the Riemannian manifold of interest. 
This can be seen from the eigenfunction expansion, however it may hold even when an eigen- 
function expansion does not exist. We assume this to be the case in the following discussion. In 
general, when we talk about the mass of the particle what we mean is \/m^ + Xmin- Since the 
Laplacian is a positive non-definite operator Xmin > on L^(A^) integrable functions, in order 
for the u-integral thus the matrix elements to be finite, energy eigenvalues are restricted to the 
following range 

-m< E < +m (49) 

The upper bound is trivially satisfied for an attractive interaction as is the case here. The lower 
bound, though, is nontrivial and indicates when our approximation of neglecting pair creation 
can no longer be trusted. If the interaction is strong enough to lower the energy by 2m from its 
non-interacting value of Hq + m then pairs will inevitably form. This is because the minimum 
energy needed to create a pair of bosons out of the vacuum is 2m. One would expect this 
approximation to break down earlier since the bosons, due to the absence of conserved quantum 
numbers, need not pop out of the vacuum in pairs. 

The lower limit of the w-integral is also finite for distinct points since the heat kernel ap- 
proaches a delta function at early times. This is precisely the reason why the on-diagonal 
elements of the principal operator matrix diverge. The well-known short time expansion of the 
heat kernel [TS| is given below 

^ oo 

lim^K^{x,x) ^ ^^^c,(a;)n'= (50) 



^We call u in Ku{x,y) the time variable, however its inverse has mass dimension 2. 
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where the coefficient functions Ck{x) are given in terms of the Riemann curvature tensor and the 
covariant derivatives evaluated at the point x, with co{x) = 1. In our case of n = 2, it is clear 
that only the first term will lead to a divergence and the divergence will be a logarithmic one. 



We see that the divergence structure is the same as what we found on in eq. (!36|) if we identify 
the inverse of the time cut-off with the momentum cut-off 1/e = A and take the inverse of the 
upper limit at which the small time approximation breaks down to be roughly ~ m — -E. 
The correspondence between the divergence structures is to be expected as renormalization 
modifies the small distance behavior of the theory and that is exactly when every Riemannian 
manifold looks like flat space. Addition of the first term in the brackets in the first line of eq. (H7j) 
removes this divergence. 

We utilized the orthofermion algebra technique to transform the original problem into one 
where we look for eigenstates of the principal operator matrix with vanishing eigenvalue. It was 
by no means necessary to use this technique, in fact we will now find the explicit form of the 
single particle state wavefunctions using a direct approach. 

3.3 Single Particle State Wavefunctions 

The 1-particle field satisfies the following operator equation 



with H given by the e — O"*" limit of eq. ([T]). The free part of the Hamiltonian suggests that we 
look for a solution of the form 



Plugging the 1-particle field ansatz above into the Hamiltonian of eq. ([I]) and grouping terms 
we find that the momentum space wavefunction must satisfy the following equation 




(51) 



(52) 




(53) 




N 



where 




(54) 
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Here, C is a normalization constant which will be determined later. Expressing V'(fli) in terms of 
this momentum space wavefunction leads to the following consistency relation among the values 
of the wavefunctioij^ at the positions of the Dirac-deltas 



1 


f \fpia^)\' ' 


9i 


Jp ^pi^p - E) 



N 



Qj f fpi(^j)fpi(^i 



(55) 



This set of equations relate the values of the wavefunction at the position of one Dirac-delta to 
the others. This equation can be written in matrix form as follows: 



<^ijip{aj) = 



(56) 



This is precisely the matrix form of the principal operator appropriate for a single boson that we 
encountered in eq. This equation does not fix the values of the wavefunction at the position 
of all Dirac-delta centers, rather it determines all the ratios. Therefore, we prefer instead to use 
the following ratios when writing down the expression for the wavefunction 



R{aj 



(57) 



As we will see, the wavefunction diverges at the positions of the Dirac-delta centers, however the 
above ratios are finite. In order for this set of equations to have solutions, the determinant of 
the principal operator matrix must be zero. This will be the case only for special values of the 
binding energy E. There will generically be different values of the binding energy for which 
this is the case. 

Normalization of the single particle state wavefunctions determines the constant C in the 
momentum space wavefunction. 



< 



> 



apilj*{p) / al^{q) 

Jq 

i2 

C y9i9j^ {aimaj) / — ^^—^ = 1 

^ JpUJ^iuJr,- E) 



I J Jp^p V^P 

When written in terms of the heat kernel the above expression yields 



(5^ 



C 



1 /"OO noo 

V ^ ti^ -^0 Jo 



a,; 



(59) 



^Even though we loosely call this the wavefunction, it should be kept in mind that it is not the generalized 
Fourier conjugate of the momentum space wavefunction because of the factor y/cJ^ in the denominator in the 
second line of eq. ([51)1 . 
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Incidentally, this equation shows that if the same problem were solved in 3 or more dimen- 
sions, a wavefunction renormalization would also be necessary. The ultraviolet behavior is worse 
than what Altunkaynak, Erman and one of us [Ij found for the non-relativistic case where the 
wavefunction normalization is necessary only in dimensions 4 and higher. When we write the 
wavefunction in configuration space we find 



N 



Using the fact that 



i=l 



hm M!l = 1 



fp{ai)fjAx) 
ujp{ijjp - E) 



(60) 



,^0+ gj{e) 

the wavefunction can be written in terms of the heat kernel as follows: 

N 



(61) 



n 



X 



-1/4 



1=1 



ds e-'^l^ / du e-^'^'e^^v^K, 



/•OO /*oo 

j,k ^0 Jo 



-1/2 



(62) 



It is worthwhile to see what the expectation values of the relativistic and interaction energies 
are for these 1-particle states. For simplicity, let's calculate the relativistic energy carried by 
these 1-particle fields for the single Dirac-delta function center case. 



<^\yj apij\p)j i^j ^,V^(g)a; ) |0 > 
Jp [Up - E) 

\fM)\" 



[^p - EY \Jq Uq [Uq - E)" 



(63) 



With the use of the heat kernel one can convert the expression for the relativistic energy into 
the form below: 



ds s^e~''/^ / du e-"™'e'^^ir„l 



ai,ai 



/"OO poo 

/ ds e-''/^ / duy/^ e-"'"'e^^^irj 
Jo Jo 



1 -1 



(64) 



The relativistic energy carried by the field is infinite! Let's also calculate the energy the field 
has due to the interaction. 



-g'^{fii)ij*{ai)ij{ai 
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LUp [ujp - E) Jg ujg [ujg - E) 

= -c-y(/.oi^(aor /-tf^ 

Jp ^p K^p - ^) 

Jp \{ujp-E) ujp[ujp-E) J 

= - < ^|i7o^ > +E (65) 

In going from the second line of this equation to the third, the relation in eq. (|55l) appropriate 
for the case of a single Dirac-delta center was used. Similarly, one obtains the last line from the 
previous one by substituting the value of the normalization constant in the single Dirac-delta 
center case. The interaction energy is also infinite and delicately cancels the infinity in the 
relativistic energy. A state with infinite relativistic energy is not actually in the domain of the 
non-interacting part of the Hamiltonian. However, the self-adjoint extension of the Hamiltonian 
extends the domain of the Hamiltonian to include such eigenstates, which would otherwise 
be excluded, with finite (total) energy. Conclusions regarding the relativistic and interaction 
energies carry over to the many center case with no modification at all but the algebra is more 
complicated. 



3.4 Bound State Energies in the Tunnelling Regime 

In the aforementioned limit of md^*" ^ 1, which hereafter will be called the tunnelling regime, 
it is possible to evaluate the off-diagonal elements of the principal operator matrix for single 
particle states analytically up to corrections of 0{{mdij)~^). Of course, on Riemannian manifolds 
dij = d{ai, aj) should be interpreted as the geodesic distance between points Oj and aj and d^*"' 
is the smallest of these geodesic distances between different points. 

-| /»00 /"OO 

= ^ e-''/^ds e-''"''Kuiai,aj)e^^'du (66) 

Jo Jo 

The equation above giving the off-diagonal elements on Riemannian manifolds provides our 
starting point. We use the following identity 

This identity expresses the fact that the heat kernel is invariant under a simultaneous scaling 
of the distances and the time variable. However, the Ihs of the equation has mass dimension n 
in n-dimensions whereas the heat kernel on the rhs is dimensionless. That's the reason for the 
prefactor on the rhs in the above equation. It is possible to obtain this prefactor by using the 
stochastic completeness condition of the heat kernel given in eq. ( l22i) . We also need the small 
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time approximation for the heat kernel on compact manifolds [20] given by 

(68) 



g-d2(x,y)/(4t) 

lim Kt{x,y) - -j—^p{d{x,y)) 



k 

Lj=o 



Moreover, we need to keep only the first term in the sum since the other terms are higher order 
in an expansion in powers of a~^. We arrive at the following results: 



p(l)co(ai,a,;(7rf-2) ^-md^.^T^iHU^^ 
$,,(E = /.,) = — ,^ X 

mdij a/1 - ifii/my 



1 {mdij)~^/'^ <^ fJ-i/m ^ 1 — {mdijY^ 



(69) 



There is an easy way to understand the parametric dependence of the exponential in the final 
result eq. (1691) . 



If {mdij)^^^'^ <^ —Hi/m < 1 the result is, except for the exponential, parametrically differential. 

^ /^(l)co(a,,a,;^dr.2) e""^'^'^ , 

$ij(E = = — — — {mdij) ' < -in m < 1 (71) 

Itt mdij\fii/m\ 

The equations given above are also valid on flat surfaces with the substitution p(l)co(aj, a^; gd~j^) = 
1. Details of this calculation can be found in the appendix. We can now use these expressions 
in eq. ( H3i) to find the change in the binding energy in the tunnelling regime due to the presence 
of the other Dirac-delta centers. 

In the tunnelling regime, we know what the ground and excited states are, these are precisely 
the eigenvectors obtained by perturbing the eigenvectors of the diagonal part of the principal 
operator matrix. We also know that these states will have energy no greater than m provided 
that fii < m and that the difference between the two is sufficiently large. This is because the 
change in the binding energy due to the other Dirac-delta centers is infinitesimally small in this 
regime and can not alter this inequality. However, in other regions of the parameter space, we 
have so far taken the existence of the ground and excited states with energy less than m for 
granted. We will now partially bridge this gap by proving that there always exists a ground 
state with energy no greater than m. 

First, we begin by showing an important result for the variation of the eigenvalues with 
energy. Eigenvalues are given by the following expectation value 

X^{E) =< x"'{EmE)x'^{E) > (72) 



■^As bi approaches zero from above and below there are transition regions where the coefficient interpolates 
between 1 and 1/2 in one case and where the parametric behavior changes smoothly in the other, so that there 
is no issue of non-analiticity. 
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The states x^{E) are eigenvectors of the principal operator with eigenvalue \"^{E), not neces- 
sarily zero. At the bound state energies, (only) the eigenvalue(s) corresponding to the bound 
state eigenvector(s) becomes zero. According to the Feynman-Helman theorem, the change in 
these eigenvalues due to a change in energy is given by the following equation 

^ r(^)I^X-(^) > (73) 

There is no summation over m in the above equations. Inserting the principal operator in eq. (1471) 
into the above equation we arrive at the following: 



dX"^ -1 



dE xAF 



/OO poo 



/*oo noo n 

Jo Jo J 



2 



-1 

< (74) 



In going from the first line to the second we have used the reproducing and symmetry properties 
of the heat kernel as well as the fact that the heat kernel is real. As we will see, what is important 
for our problem is not the sign of the derivative but that it has the same sign for all eigenvalues. 

According to the Cauchy interlacing theorem [21], if the last row and column of an (A^ + 
1) X (A^ + 1) matrix is deleted, the eigenvalues of the resulting matrix X'{E) are interlaced by 
the eigenvalues of the original one X{E) in the following manner: 

X,{E) < X[{E) < X2{E) < X',{E) < ■ ■< X'j,{E) < Ajv+i(^) (75) 

If we begin with a single Dirac-delta center with binding energy fi[ < m then adding a second 
Dirac-delta center, according to the Cauchy interlacing theorem, produces two eigenvalues which 
satisfy the following inequality 

XM) < Kif^'i) = < XM) (76) 

If we then adjust the energy to fii to make the eigenvalue of the new matrix zero i.e. Ai(/ii) = 0, 
as a consequence of eq. (174|) this new binding energy satisfies fii < fi'i < m and is the new ground 
state energy. One can repeat the same argument to conclude that there will be a ground state 
when another Dirac-delta center is added and so on. This completes the proof that there always 
is a ground state with energy not exceeding the mass of the particle. The same argument would 
have applied had the sign in eq. ([71]) been reversed but the ground state energy would have been 
/X2 satisfying A2(/i2) = 0. The third case in which the eigenvalues do not change with energy 
corresponds to a principal operator that is identically zero and is, therefore, irrelevant. 

Perron- Frobenius theorem [22] from matrix analysis guarantees that the ground state with 
energy less than m whose existence was proved by the Cauchy interlacing theorem is the unique 
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ground state. According to this theorem, a matrix A G M"^"' whose elements are all positive 
has an eigenvalue Xpp that is real and positive, with positive left and right eigenvectors, for any 
other eigenvalue A we have Xpp > |A| and the eigenvalue Xpp is simple that is it has multiplicity 
one. 

The reality property is trivially satisfied in our case, since the principal operator matrix is 
hermitian, however it does not fit the description in the theorem above. In spite of this, it can 
be rendered positive as follows: 



Here E^, is the lower bound for the ground state energy. According to the Perron-Frobenius 
theorem, the matrix M has a unique eigenvalue that is greater than all other eigenvalues, this 
implies that the matrix ($ — has a corresponding eigenvalue which is more negative 

than all the other eigenvalues and is unique. This eigenvalue flows to —^^°'^{E^,) as the energy 
is changed. During the flow it retains its uniqueness and becomes a bound state when Xpp = 
_^max ^j^^^ since the eigenvalues were all shifted down by this value. This state is the ground 
state because all the other eigenvalues intersect — at a higher value of energy due to 
the negative sign in eq. ( ff4l) . This completes the proof that the ground state is unique. 

Moreover, the fact that the elements of the left and right eigenvectors, which in our case 
are the same, can all be chosen positive means that the ground state has no nodes just as in 
elementary quantum mechanics. This is because the elements of the eigenvectors correspond to 
R{ai) in eq. (162!) and when these are all positive the wavefunction is positive over the whole 
manifold. 

4 Lower Bounds on the Ground State Energy 

The energy of the ground state should neither be too high nor too low so as not to violate 
our approximation of neglecting pair formation due to the interaction with the source. We have 
proved, in the previous section, that there always exists a ground state with energy not exceeding 
the boson mass. In this section, we will show that it is also possible to constrain the ground state 
energy from below by adjusting the parameters in the problem like the interspacing between the 
delta-function centers, their number and the minimum binding energy to a single Dirac-delta 
center. 

In order to place lower bounds on the ground state we will use a theorem from matrix analysis. 
Gersgorin theorem [22] states that the eigenvalues of the renormalized principal operator matrix 
$ are located within the union of discs 



(77) 



N 



N 




(78) 



i=l 
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where the eigenvalues are A, = in our problem. If we do not want to have any solutions below 
a given value of energy then the above inequality should not be satisfied for any i. This leads 
to the following inequality between the minimum on-diagonal and maximum off-diagonal matrix 
elements of the principal operator matrix 

N 

\^n{E,)\ > I'^nmr'' >{N-i) i$.,(i?*)r'^^ > 5^i<i>,,(e,)i (79) 



The key point here is that once we find such an E^,^ then the absence of solutions is guaranteed 
for E < E^. This is because the absolute value of the off-diagonal elements of the heat kernel 
decrease with decreasing energy for all energies whereas the absolute value of the on-diagonal 
ones increase for E < /i™'*". The ground state Egr will then be restricted to the range Egr > E^. 

4.1 Flat Space 

We found the matrix elements of the principal operator on fiat spaces. They were given by the 
following formulas: 







$n = 






27r \m — Hi J 




-1 ds 


27r Jo Vs^ + 1 



exp 



-d^j (rnVs"^ + 1 - Es^ (80) 



In order to place an upper bound on this integral we use the following relation 

< — / e-'^'^^"'-^'^'ds 



/vr _ ^ 

~ 27rdij{m- E) ^ ^ 

Therefore, there will be no solution to the eigenvalue equation for values of energy below E^: for 
which the following bound is satisfied 

>,yf-#^i) ^ (82) 

\ Jmm ^ , .mm\ \ v / 



^ . rain \ Jmin . min\ 

III /ij v'"' h^i ) 

where /i™^™ is the smallest of the binding energies to single Dirac-delta function centers and 
W{zhi{z)) = z is the Lambert W-function, also called the omega or the product-log functioE0. 
The corresponding lower bound for the ground state energy Eg^. is 

E.>ra-^a-^r)w{ ^}l;^^ (83) 

This result shows that for a fixed number N of Dirac-delta centers it is possible by choosing the 
c?™" properly, to raise the ground state energy above a critical value so that the neglect of pair 
creation is indeed justified. 



^The limiting forms of W{z) are: W{z) ~ z + 1 for z < 1 and W{z) ~ ^ {l + i^^) for z > 1. It is defined 
only for z > 0. 
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4.2 Compact Riemannian Manifolds 

The compact Riemannian manifolds over which we will investigate the issue of a lower bound 
on the ground state energy are those that have a Ricci curvature greater than or equal to zero 
everywhere. The n-dimensional sphere provides a well known example to such manifolds. 
We need to constrain the on-diagonal elements of the principal operator matrix from below and 
those of the off-diagonal ones from above. In order to achieve this, we look for upper and lower 
bounds on the off- and on-diagonal elements of the heat kernel, respectively. Below is the first 
corollary [23] that we will need. 

Corollary 1: Let be a compact manifold. Suppose that the Ricci curvature of Ai satisfies 
RicM > -K, K>0. Then W > and x E M 

K^{x,x)<y^^+A'u--/' (84) 

where n = dim Ai, A' = A'{d, K, V{Ai)) > and V{Ai) is the volume of the manifold and d 
its diameter. 

The following corollary [24J constrains the off-diagonal elements of the heat kernel from above. 
Corollary 2: Assume that for some points x,y E Ai and \/u > 0, 

C C 

Ku{x, x) < -j^ and K^{y, y) < —— (85) 

g{u) 

where / and g are increasing positive functions on (0, +oo) satisfying the regularity condition 
given below. Then, Vm > 0, i5 > 2 and for some e > 

K,(a;, y) < , ^ exp ( (86) 
regularity condition: There are numbers B > 1 and b > 1 such that 

f{s) f{u) 

VO<s<uandC = C{b, B),e = e(6, B). 

By inspection of eqs. (1841) and (1851) . one sees that these functions can be chosen as 

with C = 1. It is obvious that these functions are positive and one can easily verify that they 
are increasing. They satisfy the regularity condition fl57|) with B = 6"/^. Therefore, we can 
apply Corollary 2 to conclude that the upper bound for the off-diagonal elements of the heat 
kernel is 

K^{x,y) < C{e) + A{e){4nu)-/'^ exp (-%^) (89) 
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In order to place a lower bound on the diagonal elements of the heat kernel we make use of the 
following theorem [19]. 

Theorem 1: is a complete Riemannian manifold with Ric > 0. We have 



Kuix,y) > (47rM)-"/2exp [-d\x,y)/iAu)] 
Vx, y E M. and u > 0. Therefore, we find the lower bound to be 



(90) 



(91) 



Plugging these expressions in the integrals for the diagonal and off-diagonal matrix elements 
of the principal operator gives the following bounds 
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m — E 
m- fii 

A(e)- 



27T 
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Aie) + 



mViM.) dm 
27r 

m(m - E)ViM) 



dijim — E) 



'D/2 



dijim — E) 



(92) 



Combining these gives us the following inequahty that needs to be satisfied for a solution to the 
eigenvalue equation not to exist 



m — E^ 
m — /i™ 



> W 



4C(e)(iV-l)/D72 
mrf^*"(m - /i™")2\/(;V/l) 



[Aie)m^ViM) + tt] 



(93) 



In obtaining the inequality given above, we have replaced m — E multiplying y4(e) by 2m since 
E > —m and this only gives a more stringent constraint than the previous one. So the ground 
state energy satisfies the following inequality 



Egr>m-{m~ /i™") 



4C(6)(iV-l)/D72 
md^™(m - ^™")2y(A^) 



[A(e)mV(A<) + tt] 



(94) 



On compact manifolds, there may be more than one distance minimizing curve connecting 
distinct points. Away from a fixed point, the set of points for which this is the case are called 
the cut loci of that point. The on- and off-diagonal bounds for the heat kernel given in this 
section are valid even if the Dirac-delta centers fall on the cut loci of one another. 



4.3 Cartan-Hadamard Manifolds 

A manifold Ai is called a Cartan-Hadamard manifold [25] if is a geodesically complete, simply 
connected, non-compact Riemannian manifold with non-positive sectional curvature everywhere. 
The n-dimensional fiat and hyperbolic spaces H" are the best known examples of Cartan- 
Hadamard manifolds. The mass of the particle on such manifolds is shifted from its value on 
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compact and flat surfaces, due to tlie geometry of the manifold, to \/rn? + \min- On Cartan- 
Hadamard manifolds, the on-diagonal elements of the heat kernel are constrained from below in 
the manner given below [25] . 

Theorem 3: Assume that the sectional curvature inside the ball -B(o, r) is bounded below 
by -K^axi^)- Then, for all u > 0, x e M and 5 > 0, 

Ku{x, x) > exp [-(Ai(A^) + S)u\ (95) 

where c = c(o, 5) > 0. The minimum eigenvalue of the Laplacian Ai(A^) is restricted to the 
following range [261 EQ] 

\{n - \fKl,, > X,iM) > \in - ifKl,^ (96) 

for a Cartan-Hadamard manifold whose sectional curvature is bounded from above by —K'^^^. 

Here is an important theorem [27] regarding Cartan-Hadamard manifolds. 

Theorem 4: Any Cartan-Hadamard manifold Ai of the dimension n admits the isoperimet- 
ric function I{v) = kv~^ , k > 0. 

Another theorem [28] that we need is below. 

Theorem 5: Assume that manifold Ai admits a non- negative continuous isoperimetric 
function J(t>) such that I{v)/v is non-increasing. Let us define the function f{u) by 

/■^("^ vdv , , 

-^V„ m '''' 

assuming the convergence of the integral in the above equation at 0. If the function f satisfies 
in addition the regularity condition flHTl) then, Vx, ?/eA^,u>0,D>2 and some e > 0, 



The isoperimetric function is non-negative and the function obtained from eq. ( I97I) 

meets all the requirements of this theorem including the regularity condition (IHTI) (with any 
S > 1 in 2 dimensions). Hence, the upper bound on the off-diagonal elements of the heat kernel 
on Cartan-Hadamard manifolds is given by 

The bounds on the on and off- diagonal elements of the heat kernel lead to the following in- 
equalities for the matrix elements of the principal operator 



2tt \mcH - f^i 
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where 

m'cH -rn' + \kI^ + 8{K') (102) 

Combining these gives us the following inequahty that needs to be satisfied for a solution to the 
eigenvalue equation not to exist 



m-E. / (iV-l)C(.,K)v/IV2 

I „/r\ Jmin / ..min\ I V / 



mcH - I^T^'^ \c{5)dlj'''''' {mcH - l^i 

We replaced uich in the numerator on the Ihs by m as mcH ^ fTi and this replacement only 
results in a more stringent constraint than the one before. Then the inequality for the ground 
state becomes 



F >m (m ^ w( - >^)Vm \ 



(104) 



Using the large z approximation of the Lambert W-function, one can check in the perturbative 
limit, that the ground state energy on fiat and compact manifolds approaches m logarithmically 
from below. However, on Cartan-Hadamard manifolds, the bound we obtain for the ground 
state in the perturbative limit is not as stringent. We could only prove that Egj. > m whereas 
on such manifolds we would expect Eg^ > y^m^ + Xi(M.) > m. 



5 Conclusion 

In this paper, we utilized the orthofermion algebra technique to convert the problem of deter- 
mining the spectrum of relativistic bosons interacting with N localized sources to one where we 
had to solve for the eigenstates of the finite dimensional principal operator matrix with vanishing 
eigenvalue. The Delta-function type Hamiltonian that we are interested in was not well-defined 
and hence would lead to a principal operator matrix that was not well-defined either. We ren- 
dered the original Hamiltonian well-defined by the use of the heat kernel which approaches a 
Dirac-delta function in the limit e — > 0"^. 

The principal operator was found using this Hamiltonian and it was noticed that in order 
to get physical results that did not depend on the cut-off, the coupling constant had to depend 
on the cut-off in a particular manner. After regularizing the principal operator, we chose the 
renormalization condition and made the natural choice for which the bound state energy in the 
case of a single Dirac-delta function center case coincided with the experimentally measured 
binding energy to that center. 

We could not find the eigenstates of the principal operator and the bound state energies 
corresponding to them in the general case. However, in the tunnelling regime, we managed to 
find the eigenstates and the bound state energies on 2 dimensional compact manifolds and flat 
surfaces. Also, we found the wavefunctions for single particle states on general manifolds and 
determined that multiparticle states could be written as products of these single particle states. 
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We made use of the Cauchy interlacing theorem to prove that a ground state with energy not 
exceeding the boson mass exists so long as the single Dirac-delta center binding energies are not 
above this mass. Most importantly, we proved that the ground state energy has a lower bound 
on flat, compact and Cartan-Hadamard type manifolds for appropriate choices of the number of 
Dirac-delta centers N, their separations dij and the initial binding energies fii to single Dirac- 
delta centers. Thus, validating the neglect of particle production due to the interaction with the 
source, on a large class of 2-dimensional Riemannian manifolds. 

We would like to extend the proof that the ground state with energy no greater than the 
boson mass exists to excited states or to show under what conditions excited states that meet 
this requirement exist in future research. It would also be interesting to find approximate/exact 
solutions in regimes where the bound state energy is fairly different from the binding energy to 
the single Dirac-delta center. 
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7 Appendix: Calculation of Bound State Energies in the 
Tunnelling Regime 

In this appendix, we give those details that were omitted in the main body of the text related 
to the calculation of bound state energies in the tunnelling regime. We first consider the case 
^-1/2 ^ jj. ^ I ^ Q-i where a = mdij and hi = /ij/m. We start with eq. (166!) . 

-| /"OO /»oo 

1 nco POO 2 

= — ^/ ds rfMe-""'i^„(ai,a,) e~M^^) e"^' 
Jo Jo 

= Hds r du K^{a,,aA e'l-^v^^)'/^ e-"'"'(^-^'), 6 = -, 

VTT Jo Jo rn 
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then we plug the following equation in the one above 

1 



4 



and we get 



^1 



2 du 1 (n n ■ n^-^\ -u(l-h'')m^/A f 



ds 



A 

oo 



z-'^ds 



where we have chosen A = -j^. At this point, we argue that the s integral is independent of u 

"■ij 

provided b^/ua ^ 1. This will be the case if the u integral is dominated by w's of 0{1). This 
turns out to be case as can be checked a posteriori. The other is a necessary condition that 
/li must satisfy in order for this method to work. Another condition that has to be satisfied in 
order for this to be a good estimate to the integral is that the exponential be highly damped 
that is a (1 — 6^) ~ a(l — 6) ^ 1 for 6 ~ 1. For values of u that dominate the integral, the time 
argument of the heat kernel is a small number and the small time approximation for the heat 
kernel given below may be used. 



g-d2(x,J/)/(4t) 

Um Ktix.y) ~ — -- — , ,^ p(d(x,y)) 



J2t'cj{x,y;g) + 0{t'+') 

.3=0 



Once the small time approximation is written down it can be seen that only the first term should 
be kept since the higher order terms are suppressed by powers of a~^. We continue with the 
following equation after having taken the above simplifications into account 
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The huge damping of the exponential allows us to make a saddle point approximation. We need 
the value of the function and its second derivative at the point where it is stationary 
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using these in the equation above we find the following 
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27r 



which is eq. (169|) if we make the identification 6 = 6j. The case ^ a~^/^ can be calculated in 
the same way except that the s integral gives a different factor. The case —b ^ a~^/^ requires 
a different approach than the previous one. We begin again with eq. (1661) 
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Since this exponential is highly damped the major contribution will be from the t ^ region 
where we can ignore the t"^ term in the exponential and in the denominator. These simplifications 
lead to the elementary integral and the result given below 
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which is eq. (I7T1) with b = bt. 
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